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ABSTRACT 

We present the results of deep, high angular resolution Very Large Array 
(VLA) surveys for radio continuum emission towards a sample of 14 GLIMPSE 
Extended Green Objects (EGOs). Identified as massive young stellar object 
(MYSO) outflow candidates based on their extended 4.5 ^m emission in Spitzer 
images, the EGOs in our survey sample are also associated with 6.7 GHz Class 
II and/or 44 GHz Class I CH3OH masers. No continuum is detected at 3.6 or 
1.3 cm towards the majority (57%) of our targets (median rms ~0.03 and 0.25 
mJy beam~^). Only two EGOs are associated with optically thin emission con- 
sistent with ultracompact /compact HII regions. Both of these sources exhibit 
cm— A multiplicity, with evidence that one of the less-evolved members may be 
driving the 4.5 /um outflow. Most of the other cm-A EGO counterparts are weak 
(< 1 mJy), unresolved, undetected at 1.3 cm, and characterized by intermedi- 
ate spectral indices consistent with hypercompact (HC) HII regions or ionized 
winds or jets. One EGO cm counterpart, likely an optically thick HC HII re- 
gion, is detected only at 1.3 cm and is associated with hot core line emission 
and H2O and 6.7 GHz CH3OH masers. The results of our exceptionally sensitive 
survey indicate that EGOs signify an early stage of massive star formation, be- 
fore photoionizing feedback from the central MYSO significantly influences the 
(proto) cluster environment. Actively driving outflows (and so, presumably, ac- 
tively accreting), the surveyed EGOs are associated with significant clump-scale 
gas reservoirs, providing sufficient material for sustained, rapid accretion. 
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1. Introduction 

A key difference between tiigli-niass stars and ttieir low-mass counterparts is tliat mas- 
sive stars are sufficiently luminous to emit ionizing radiation. Early efforts to compile 
large samples of massive young stellar objects (MYSOs)0 focused on ultracompact (UC) 
HII regions, identified by their cm-A radio continuum emission and observationally defined 
as h aving sizes < 0.1 pc , dens i ties > 10^ cm~^ 



and emission measures > 10'' pc cm ^ 



[e.g. IWood fc Churchwelll 1 19891 : iKurtz et al.l Il994j : iHoare et al.l 120071 ). The main mech- 
anisms proposed to explain why UC HII regions are not always observed around lumi- 
nous (^lO'' L0) MYSOs-" quenching" , gravitational trapping, and swelling- are consequences 



Keto 



2007 



Hoare &: Franco 2007; Yorke fc Bodenheimer 2008: Hosokawa &: Omukai 



Hosokawa et al.l I2OIOI ). This led iHoare fc Francol (120071 ) to suggest that UC HII regions 



of the high accretion rates assoc i ated with massive star formation (e.g Walmslev 19951; 

20091 : 



the transition is not clear-c ut (e.g 



Keto 


2002; 


Keto & Wood 


2006: 


Keto & Klaassen 


2008 



Galvan-Madrid et al.ll2009l ). While an evolutionary sequence for MYSOs is still under de- 
velopment, it is generally accepted that a UC HII region represents the photoionization of 
its environment by a young or early B type star, and thus a late stage of m assive star 
formation (e.g. IZinnecker fc Yorkd 120071 ; lEUingsen et al.l 120071 : iBreen et al.ll2010l ). A good 
example of this st age is G5. 89— 0.39, where the can didate ionizing star has been identified 
(IFeldt et al.ll2003l ) in the midst of a dust-free cavity (IHunter et al.ll2008l ). 



The first appearance of observable emission from ionized gas associated with a MYSO is 
considerably less well understood. Hypercompact (HC) HII regions, ionized accretion flows, 
outflows, jets, stellar winds, and ionized/photoevaporating disks all may exist around form- 



density 



> 



10^ cm-3 



(e.g. 



on," 


remain the subiect of considerable discussion (e.g. 


Menten & van der Tak 


2004 




Keto 


2007 




Gibb & Hoare 


2007 




Hoare & Franco 


2007: Keto & Klaassen 


20081: 


Rodrfeuez et al. 


2008 


Galvan-Madrid et al. 


I2OO9). HC HIIs are defined by their scale: size 5, 0.05 pc. 



Hoare et al.l 120071 : iLizanol 120081 ) . An additional distinction is 



that sources in which photoionization by the central MYSO dominates, including ionized 
accretion flows or disks, are considered HC HIIs. In contrast, in jets and winds, shocks 
may contribute significantly to t he ionization (see for example IShepherd et al.l 120041 : iKeto 
20071 : iGalvan-Madrid et al.ll2010l ). HC HIIs, winds, and jets can be difficult to distinguish 
observationally, as all are characterized by intermediate cm-A spectral indices (S^, oc z/". 



^We define MYSOs as young stellar objects (YSOs) that will become main sequence O or early-B type 
stars (M>8 Mq). 
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— 0.1< a < +2) and all are generally fainter than UC HIIs at cm wa velengths ( jSewilo et al. 
20041 : iKurtzlboOsi Icibb fc Hoar^boOTl : iKeto et alj[2008l : iLizanolboOsh . Thus the distinctions 
among these phenomena are in large part dynamical: a key observational discriminant is 
the recombination line width (or for jets, proper motion), with jets having the l argest veloc- 



ities/linewidths, followed by winds, then HC HII and finally UC HII regions (IHoare et al. 
200?! : iHoare fc Francol [20071 ) . 



Despite the uncertainties outlined above, the strong association between ongoing ac- 
cretion (and hence youth) and absence of a UC HII region means that characterizing the 
presence/absence of cm- A continuum emission is a key step in understanding the evolu- 
tionary state of any sample of MYSOs. A large new sample of MYSOs with active out- 
flows-and so presumably ongoing accretion-has been identified from Spitzer Space Telescope 
surveys of the Galactic Plane based on extended 4.5 a m emission. Kno wn as "Extended 



Gree n Objects (EGOs)" (ICvganowski et al.ll2008l . l2009l ) or "green fuzzies" (IChambers et al. 



20091 ) for the common coding of three-color IRAC images, their extended 4.5 fim emission 



is thought to arise from molecu 


ar line 


emission shock-excited in prot 


□stel 


ar outflows (pre- 


dominantly Ho: 


Smith & Rosen 


2005 




Smith et al. 


2006; 


Davis et al. 


2007 




Ybarra & Lada 


2009; 


Ybarra et al. 


2010; 


De Buizer & Vacca 


2010|). 


Cyganowski et al. 


(2008 


, hereafter COS) 



cataloged over 300 EG Os in the Galactic Lega cy Infrafred Mid-Plane Survey Extraordinaire 
(GLIMPSE-I) survey (IChurchwell et al.ll2009l ). and argued that their MIR properties were 
consistent with young, embedded MYSOs. To tes t the hypothesis that GL IMPSE EGOs 
are specifically massive YSOs with active outflows, ICyganowski et al.l (120091 . hereafter C09) 
conducted a VLA search for two diagnostic types of CH3OH masers t owards EGOs: 6.7 



GHz Class II CH3OH masers, associated exclusively w ith MYSOs (e.g. iMinier et al.l 12003 



Bourke et al.ll2005l : IXu et al.ll2008l : iPandian et al.ll2008h . and 44 GHz Cl ass I CH.sOH masers 



assoc i ated with molecul ar outflows and outflow-cloud interactions (e.g. iPlambeck fc Menten 



I990I : iKurtz et al.l 120041 ). The detection rates were remarkably high, confirming the MYSO 
and outflow nature of the EGO targets: of a sample of 28 EGOs, > 64% have associated 6.7 
GHz CH3OH masers, and ~89% of the 6.7 GHz maser sources also have associated 44 GHz 
masers. A c ompl ementary JCMT survey f ound evidence for warm dense gas associated with 
EGOs (IC09l ) and I Cyganowski et al.l (120 111 , hereafter Cll) found hot cores and high- velocity, 
bipolar molecular outflows in two EGOs studied in detail with high-resolution mm observa- 
tions. All of this is strong evidence for EGOs being young, actively accreting MYSOs, but 
constraints on the cm-A continu um p roperties of EGOs are conspicuously lacking. The 44 
GHz continuum observations of |C09| ruled out bright UC HII regions as powering sources 
for most (95% of) targeted EGOs, but the sensitivity was insufficient to detect fainter UC 
HII or HC HII regions. In this paper, we pres ent d eep, high-resolution VLA 3.6 and 1.3 cm 



continuum observations of 14 EGOs from the |C09| sample, with the aim of constraining the 
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evolutionary state of EGOs. 

2. Very Large Array Observations and Data Reduction 

We observed a sample of 14 EGOs with the Very Large Array (vlaH at 3.6 cm (8.46 
GHz) and L3 cm (22.46 GHz) using the standard continuum mode (2 x 50 MHz, dual 
polarization, for an effective bandwidth of 172 MH^). The 3.6 cm data were taken in the 
VLA B configuration on 2009 May 7 and 14. Each day, short observations of each target were 
staggered throughout the observing block to optimize the uv-coverage, for a total on-source 
time of ~45 minutes per source. Calibration followed standard procedures in AIPS, including 
the use of a model for 3C286, the absolute flux calibrator. The quasars J1751+096, J1733- 
130, J1925+211, and J1911-201 were used as phase calibrators and to correct for closure 
errors on EVLA-VLA baselines using the AIPS task BLCAL. We estimate that the absolute 
flux calibration is accurate to ~5%, and that the absolute positional uncertainty is ~0'.'l. 

We chose the VLA C configuration for our 1.3 cm observations in order to approximate 
as closely as possible the uv coverage and synthesized beamsize (~1") of the B configuration 
3.6 cm data. The 1.3 cm data were taken on 2009 July 10 and 17 using fast switching (cycle 
time 2.5 minutes) and reference pointing. The 22 GHz zenith opacity was ~0. 16-0. 17 on July 
10 and ~0. 13-0. 14 on July 17. Cycles on each source were distributed across the observations 
to improve uv-coverage; the total on-source integration time was ~25-30 minutes per target. 
Table [Ulists the fast-switch phase calibrators. Observations of the bright quasars J1751-I-096 
and J2253-t-161 were used to correct for EVLA-VLA closure errors using the AIPS task 
BLCAL. The data were calibrated in AIPS following standard high-frequency procedures, 
including using a model for 3C286. For the July 17 data, 3C286 was used for absolute flux 
calibration. On July 10, the 3C286 data were unusable, so the absolute flux scale was set 
using J2253+161, assuming S(22.46 GHz)=6.22 Jy based on VLA monitoringS Comparison 
of the derived fluxes for the phase calibrators and J175H-096 between our datasets and with 
VLA flux monitoring shows no systematic offset, and suggests the absolute flux calibration 
is accurate to ~10%. 

The data were imaged in AIPS, and, for sources with sufficient signal-to-noise, self- 



^The National Radio Astronomy Observatory operates the VLA and is a facility of the National Science 
Foundation operated under agreement by the Associated Universities, Inc. 

■^httpi/Zwww. via. nrao.edu/astro/guides/exposure/ 

^http:/ /www. aoc.nrao.edu/~smyers/calibration/2009/K_band_2009.shtml 
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calibrated. For each target, we also convolved the 3.6 cm image to the 1.3 cm resolution, to 
facilitate direct comparison of fluxes and limits between the two datasets. Image parameters 
for each EGO are presented in Table [TJ Some fields are affected by image artifacts from 
poorly-sampled large-scale emission. In these cases, images were made first with all data to 
search for extended emission ( §4.3p . and then with uv limits to reduce artifacts and improve 
sensitivity to compact emission; uv weighting adjustments are indicated in Table [1] The 
largest angular scale of the observations is ~20", but we cannot recover fluxes for objects 
more extended than ~10". The primary beam (FWHM) of the VLA is 5.3' at 3.6 cm and 
2.0' at 1.3 cm. All measurements were made from images corrected for the primary beam 
response. 



3. Results 

Of the 14 EGOs in our sample, 6 are detected in cm continuum emission. We de- 
fine a cm-A EGO "detection" as >4a emission within the polygonal aperture for the EGO 



pubhshed by ICOSl . Table [2] lists observed properties of cm continuum sources; for EGOs 
associated with multiple cm sources, they are designated -CMl, -CM2, etc., in order of de- 
scending peak intensity. The vast majority of detected cm sources are unresolved at the ~1" 
resolution of our observations. For sources detected with sufficient signal-to-noise, a single 
two-dimensional Gaussian was fit, and the fitted centroid position, peak intensity, integrated 
flux density, and deconvolved source size are reported in Table [2l For weaker sources for 
which the size is not well-constrained. Table [2] lists the fitted centroid position and peak 
intensity from a Gaussian fit assuming an unresolved point source. For marginal detections 
(< 5a), the position and intensity of the peak pixel are reported in Tabled and the quoted 
uncertainty is the la rms in the image. The parameters of the one resolved, irregular cm 
source (G49.27— 0.34-CMl) were mea sured using the pol ygon aperture photometry program 



developed for the CORNISH survey (jPurcell et al.ll2008f jl and Table [2] reports the position 



of the brightness-weighted source center, the peak intensity, the integrated flux density, and 
the equivalent source diameter. Table [3] lists 4a" upper limits for EGOs not detected in cm 
continuum emission. Figure [1] shows the spatial relationship of the MIR, cm-A continuum, 
and Class I and Class II CII3OH maser emission for EGOs with cm-A detections. For all 
sources except Gil. 92— 0.61 (which is detected only at 1.3 cm), contours of the 3.6 cm im- 
age are shown in Figure [H Most EGO cm counterparts are detected only at 3.6 cm; for 
sources that are detected at both wavelengths, the 1.3 cm images (not shown) have similar 
morphology, but lower signal-to-noise. 



■''http://www.ast.leeds.ac.uk/~phycrp/software.html 
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In addition to the EGO counterparts listed in Table |2l we detect other continuum 
sources within the VLA field of view that are not "associated" with the target EGO (as 
defined above). Such sources are listed in Table H] if they are within a 1' radius of the 
VLA pointing center (corresponding to the FWHM of the VLA primary beam at 1.3 cm), 
so that a 1.3 cm fiux or limit can be measured. Figure [2] shows three-color IRAC images 
of a 2' square field, centered on the VLA pointing center, for each of the EGO fields i n 
Table H overlaid with contours of the VLA 3.6 cm, MIPSGAL 24 urn ficarev et allboogh . 



and Bolocam Galactic Plane Survey (BGPS; Aguirre et al. 2011 : Rosolowsky et al. 201ol ) 1.1 



mm continuum emission. The information in Table H] is equivalent to that in Table [2|, and, 
for compact sources, measured as described above. For extended sources poorly imaged by 
the interferometer, only estimated positions, sizes, and angular separations from the target 
EGO are listed in Table |H The sources in Tableware designated "F" (for "field"), the 
name of the targeted EGO, and -CMl, -CM2, etc. in order of descending peak intensity. 
The angular separations of these sources from the target EGO positions range from 8-51". 
Some are likely associated with the same star-forming region as the EGO, and are noted in 
Section [3l2l In general, however, insufficient information is available to determine association 
and/or the physical properties of these sources, and we do not include them in our analysis. 



3.1. Physical Properties: Estimates and Limits from the Cm- A Continuum 



For EGOs undetected at both 3.6 and 1.3 cm, we adopt two approaches to constraining 
the physical properties of the (outfiow) driving source and of any ionized gas. The number 
of ion izing photons corresponding to an unresolved 4a source is given (following iKurtz et al. 
mi by: 

1 \/ ^°'^M(^J^y 
Jy) \kpc) 



N'c < 7.59 X 10^® T7°-^ 



(1) 



a(z/, Tg) J \GHz 

where A*"^ has units of (photons) s~^, Te is th e electron temperature in K, a(z^, Te) is a factor 
of order unity, interpolated from the tables in lMezger fc HendersonI ( 119671 ). and units for the 
other quantities are as indicated. This calculation provides an upper limit to the ionizing 
radiation that could be emitted by the MYSO and yet be undetected in our observations. 
However, the upper limit obtained from equation [1] may be an underestimate, as it assumes 
optically thin, ionization-bounded free -free emission an d no absorption of ionizing radiation 
by dust (see for example discussion in lKeto et al.ll2008l ). 



If, in contrast, the emission is assumed to be optically thick, an upper limit to the source 
size can be derived from Si, = B^{Te)fls where S,^ is the observed limit on the fiux density. 
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B(z/, Te) is the Planck function, and Qg is the source sohd angle. In the Raleigh- Jeans limit, 



iAU] < 6.639 X 10^ 



0.5 



(2) 



where r is the radius (assuming spherical geometry), 8,^ is the flux density (limit) in Jy, D is 
the distance in kpc, Tg is the electron temperature in K, and u is the observing frequency in 
GHz. Recent studies assume electron temperatures in the range 8000-10000 K for H C and 
UC HII regions (e.g. IPandian et aPboiol : iBattersbv et al.lboiol : Ijohnston et aPboogh . The 
difference in the values of the derived parameters is <15%. We adopt Te=10000 K here. 

For all EGOs undetected at both 3.6 and 1.3 cm in our observations. Table [5] lists the 
limits on the source size and ionizing photon rate. The 3.6 cm data are significantly more 
sensitive, with a median 4(7 limit on Niyc of 2.24 xlO^"^ s~^, compared to 1.81x10"^^ s~^ for 
the 1.3 cm d ata. We incl ude both estimates because the turnover frequency is a function of 



density (e.g. iKurtzl |2005[ ) . so emission from dense HII regions is more likely to be optically 



thin at shorter wavel engths. However, emission from HC HII regions may be optically thick 
even at 1.3 cm, and iKeto et al.l ( l2008l ) discuss ways in which density gradients within HII 
regions may "stretch" the transition region of the spectrum. 

For EGOs detected at either 3.6 or 1.3 cm. Table E] lists the 3.6 to 1.3 cm spectral 
index or limit, and the estimates for the ionizing photon flux and source size calculated from 
equations [1] and |2] and the flux densities in Table [2l Spectral index limits are calculated using 
the 4cr upper limit at the undetected wavelength. The size estimates from equation |2l which 
assume optically thick emission, are omitted in cases where the spectral index indicates the 
emission is optically thin [a ~ —0.1), but retained for sources with intermediate spectral 
indices for co mparis on. For all sources, physical properties were calculated assuming the 
distances from ICOQI . which are listed in Tables [5] a nd [61 Throughout §3.2[ spectral types 
corre sponding to N f,,,^ are taken from ISmith et al.l (j2002[ ) for stars earlier than B1.5V, and 
from iDoyoru (119901 ) for later B-type stars. 



3.2. Notes on Individual Sources 



3.2.1. GlO.29-0.13 



This EGO is unique among the IC09l sample in being associated with both 6.7 GHz Class 
II and 44 GHz Class I CH3OH masers, but lacking a discrete MIPS 24 yum counterpart. The 
detectability of a weak 24 yum counterpart is, however, limited by the point-spread function 
(PSF) wings of the adjacent MIR-bright complex. This saturated MIPS 24 /im source. 
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located east of the EGO, is detected at both 3.6 cm and 1.3 cm in our VLA images (Table IH 
Fig. |2]). However, the cm- A emission is extended, and not well-imaged. Lower-resolution 
(37x25") 20 cm VLA ima ges show that this c m source is part of an extended ionized complex, 
designated G10.30— 0.15 ( iKim fc Kodl200ll ). The velocity of th e molecular gas associated 
with this ionized complex (v2,5i^=13.5 km s~^. lKim fc Kodl2003l ) agre es well with the radi o 
recombination line velocity reported for GIO. 315— 0.150 (12 km s~^, iDownes et al.lll980l ). 
and with the vlsr of the dense gas associated with the EGO (13.6 km s~^. ICOOl ). indicating 
that the EGO is part of the same star-forming region. 



3.2.2. Gil. 92- 0.61 



The EGO Gl 1.92—0.61 is the only source in our VLA continuum survey to be detected 
only at 1.3 cm. The 1.3 cm detection of 1.07 mJy beam~^ is marginal (~4.2 a), but intrigu- 
ing because it implies a spectral index >1.7. As shown in Figure [H the 1.3 cm continuum 
emission is coincident with the 1.4 mm dust continuum emission from the hot core MMl. 
The 1.3 cm peak i s ^0"4 (> 1500 AU) south west of the CARMA position of MMl Jcnl) 



and a H2O maser (IHofner fc Churchwelllll996l ). which are coincident within astrometric un- 
certainties. The cm peak is ~0'.'2 (> 930 AU) northwest of the intensity- weighted 6.7 GHz 
Class H CH3OH maser position from |C09| (Fig. [T]). The 1.3 cm detection is too strong to be 
pure dust emission, as this would account for at most 0.2 mJy (extrapolating from the com- 
parable resolution CARMA data with a spectral index of 3). If the 1.3 cm emission is due to 
optically thick free-free emission (a=2) and remains optically thick into the mm, then free- 
free emission could potentially account for ~59% of the CARMA 1.4 mm flux of MMl, with 
the remainder due to dust emission. The calculated ionizing photon flux, ^Lyc >1.66xl0^^ 
s~^, would correspond to a single ionizing ZAMS star of spectral type B2. However, there 
is strong evidence that G11.92— 0.61-MMl is actively driving an outflow, including Class 
I CH3OH masers and SiO, 4.5 /im, and high-velocity ^^C0(2-l) and HCO"'"(1-0) emission 
(jCllI : IC09l )-and therefore actively accreting. It is therefore very unlikely that the driving 
source is in fact in a ZAMS configuration fe.g. iHoare &: Francdl2007l : lYorke &: Bodenheimei 
20081 : iHosokawa fc Omukail[2009l : iHosokawa et al.lboioh . 



3.2.3. G 18. 67 +0.03 

The EGO G18.67+0.03 is the easternmost of four MIPS 24 /im sources associated with 
a ridge of 1.1 mm emission (Fig. [2]). The three eastern 24 /im sources are all associated 
with 44 GHz Class I CH3OH masers; two of these sources, including the EGO, are also 
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associated with 6.7 GHz Class II CH3OH maser s. (T he westernmost 24 /iin source falls just 
outside the field searched for maser emission by IC09l .) Interestingly, the only cm continuum 
emission in the field is detected towards the 24 /im source associated with Class I but not 
Class II CH3OH masers. Maser evolutionary sequences have generally placed Class I CH3OH 
masers among the earliest observable signs of massive star formation , due largely to their 
association with m olecular outfiows and outfiow/cloud interactions (e.g. lEUingsen et al.ll2007l : 



Breen et al.l 120101 ). Indeed, one of the brightest known Class I maser sources is NGC 6334 



I(N), which cont ains a rich cluster of c ompact millimeter continuum sources (IBrogan et al. 



20091 ) . However, IVoronkov et al.l ( I2OIOI ) have recently suggested that Class I CH3OH masers 



may also be excited by expanding HII regions driving shocks into the surrounding molecular 
cloud, and so (also) trace a much later stage of massive star formation. Because most 
high-resolution Class I maser searches to date have targeted Class II masers, there are few 
Class I CH3OH maser so urces associated with cm- A continuum emission and known to lack 
Class II CH3OH masers (IVoronkov et al.ll2010l ). The cm continuum source F G18.67+0.03- 
CMl is an important addition to this sample. Since the focus of this paper is on the cm-A 
properties of EGOs, we defer further analysis of F G18. 67-1-0. 03-CMl to a future paper that 
will present high resolution mm-A molecular line and continuum data, which are required to 
better characterize the physical properties and relative evolutionary states of the EGO and 
the other 24 //m sources. 



3.2.4. G 19. 36- 0.03 

Two weak (< 1 mJy) 3.6 cm sources are detected in this field, neither coincident with 
the EGO. The stronger cm source, F G19.36-0.03-CM1, is ~8'.'0 (> 19400 AU) northwest 



of the nominal EGO position cataloged by ICOSl . arid_~4'.'5 (>10900 AU) northwest of the 



intensity- weighted 6.7 GHz CH3OH maser position (IC09l ). This cm source is coincident with 
multiband IRAC and with MIPS 24 /im emission; a short line of 44 GHz Class I CH3OH 
masers originates near the cm source and extends ~2" NNW. F G19. 36— 0. 03-CMl is thus 
another candidate for Class I CH3OH masers possibly excited by an expanding HII region. 
The NW-SE extension of the MIPS 24 fim emission indicates that it is likely a blend of 
emission from multiple sources, possibly F G19. 36— 0. 03-CMl and the MYSO pinpointed 
by the 6.7 GHz CH3OH maser. These sources may be two members of a (proto) cluster, 
possibly in different evolutionary stages. We defer further discussion of F G19. 36— 0. 03- 
CMl to a future paper that will present high-resolution mm-A line and continuum data to 
better constrain the properties of the (proto) cluster members. The nature of the weaker cm 
source, F G19.36— 0.03-CM2, and its relation (or lack thereof) to the EGO are unclear. F 
G19.36— 0.03-CM2 is coincident with compact, multiband IRAC emission, but lacks strong 



24 fim emission (Fig. [2]), and could be a much more evolved object. 



3.2.5. G 24. 94+ 0.07 

The 3.6 cm EGO counterpart EGO G24.94+0.07-CM1 is coincident with both 6.7 GHz 
Class II CH3OH maser and MIPS 24 /im emission (Fig. H]). The fitted position of the cm 
source is offset by only ~0'.'4 (> 1100 AU) from the intensity- weighte d 6.7 GHz maser position 



and by ~1'.'4 (> 3400 AU) from the nominal MIPS 24 /im position ( lC09l ). The nondetection 
at 1.3 cm implies a <0.7. The 3.6 cm emission is unresolved, and the ionizing radiation flux 
is low (4.51 xlO"^^ s-\ ZAMS star of type B2 to B3). Since the central MYSO is actively 
driving an outflow-as indicated by the 4.5 /im emission, Class I CH3OH masers, and SiO 
emission-and accreting, however, it is unlikely to have contracted to a ZAMS conflguration 
(see also §3X2]). 



3.2.6. G25.27-0.43 

Two weak (< 1 mJy) 3.6 cm sources are detected in this fleld, both >0.5 pc in projection 
from the EGO. The weaker, F G25.27— 0.43-CM2, is coincident with the center of a nebula 
that is bright in all IRAC bands and at 24 fim (Fig. |2]), and corresponds to the IRAS source 
18362-0703 (~28" south of the EGO). The brighter cm source, F G25.27-0.43-CM1, does 
not appear to have a MIR counterpart, and its nature is unclear. 



3.2.7. G 28. 28- 0.36 



A weak 3.6 cm source, EGO G28.28-0.36-CM1, is marginally detected (~4.1a) ~1'.'2 
(> 3800 AU) northeast o f the 6.7 GHz CH3OH maser associated with the EGO (intensity- 
weighted maser position, ICOQI ). The limit on the spectral index (<1.6) is unconstraining. 



Three other cm sources are detected in the fleld. F EGO G28.28— 0.36-CMl is the well- 



Walsh et al.lll998: 


Lonemore et al. 


2007; 


Churchwell et al. 



2OIOI ). IChurchwell et al.l (120 lOl ) report a H30a recombination line velocity of 40.8 km s ^ 
(FWHM 31.5 km s'^) and a HCO+(3-2) velocity of 47.1 km s'^ towards the UC HII re- 
gion. The latter is i n goo d agreement with the molecular gas velocity measured towards the 



EGO (49.5 km s-MC09h. The EGO and the UC HII region (> 0.5 pc away) are thus likely 
part of the same star-forming complex, which may also include the compact cm source F 
G28.28-0.36-CM3 (coincident with a MIR nebula ~24"> 0.4 pc southeast of the UC HII 



region). The other field source, F EGO G28.28— 0.36-CM2, does not have a MIR counterpart. 



3.2.8. G 28. 83- 0.25 



This EGO hes at the edge of the MIR bubble N49 fIChurchwell et al.ll2nn6[ ). and several 
auth ors have investigated the possibility of triggered star formation around the bubble rim 



(e.g. iWatson et al.ll2008t IZavagno et al.ll2010l : iDeharveng et al.ll2010l ). The 3.6 cm emission 
~13" (> 0.3 pc) north of the EGO (F G28.83-0.25-CM1, over-reso lved in our VL A i mage) 
corresponds to a compact HII region seen in MAGPIS 20 cm data. IDeharveng et al.l (120 10[ ) 
quote a 20 cm fiux d ensity < 0.1 Jy and a spect ral type for the ionizing star of BOV. 
The CO (3-2) maps of iBeaumont fc Williams! ( l2010l ) show a velocity gradient of only a few 
km s~^ in the molecular gas around the bubble rim, confirming that the three 24 sources 
in F igure [2] (the compact HII region, the EGO, and the 24 source south of the EGO, see 
also ICOQI ) are part of the same star-forming clump. 



Two faint (< 0.5 mJy) 3.6 cm continuum sources are detected coincident with the 
EGO. The nature of the brighter source, EGO G28.83— 0.25-CMl, is somewhat pu zzling : 
located ~5'.'2 (>26200 AU) west of the 6.7 GHz Class II CH3OH maser position Jcooh . 
CMl is similarly offset from the 24 /xm peak and appears to coincide with a deficit of 4.5 
fim emission (Fig. [T]). G28.83— 0.25-CMl has an intermediate spectral index (<1.1) and 
modest ionizing photon flux (7.95x10^^ s^^)- The fainter 3.6 cm continuum source, EGO 
G28.83-0.25-CM2 (~4.7a), is coincident with the 6.7 GHz CH3OH maser and with 24 /im 
emission. The spectral index is poorly constrained (<1.9), and the size implied in the case 
of optically thick emission is small (~90 AU). 



3.2.9. G35.03+0.35 



The EGO G35.03+0.35 is remarkable among our sample for its cm-A multiplicity: flve 
compact 3.6 cm sources are d etected, with projected separations <20000 AU (Fig. [1]). In 
addition, iBrogan et al. I ( 1201 ll ) flnd evidence for an additional source NW of CMl and CM2, 
detected in NH3 and 1.3 mm continuum emission, but not at cm wavelengths. CMl, CM2, 
and CM3 all lie along a ridge of 24 ^m emission. T he M IPS image is saturated, so the 
position of the 24 ^um peak is highly uncertain (see also lC09l ): however, the peak of CM2 lies 
within the saturated MIPS pixels, while the peaks of the other cm sources do not. CM2 i s 
also coincident with blueshifted 6.7 GHz Cl ass II CH3OH mas ers (Av ~6-12 km s~^,|C09|), 
blueshifted OH masers (Av ~3-14 km s~^. lArgon et al.ll2000[ ). redshifted H2O masers (Av 
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13-17 km s ^. iForster &: Cas'welllllQQQl ). and hot molecular gas seen in NH3(6,6) emission 



(Tfc=220 K. lBrogan et al. 1 1201 If ) 



CMl was detected in the iKurtz et al.l (119941 ) survey of UC HII regions, and also at 
7 mm by |C09| . T o bet ter constrain the SED of CMl, we reimaged the 7 mm continuum 
data presented by IC09l . applying a uvtaper to approximate as closely as possible the beam 
of our 3.6 and 1.3 cm data. We then convolved the tapered 7 mm image with the beam 
listed in Table [H and corrected for the more recent VLA flux scale used for the 3.6 and 1.3 
cm data. The flux density of CMl measured from this new 7 mm image is 10.9±3.0 mJy 
(including a 20% uncertainty in the absolute flux calibration). In addition, we measured 
the 6 cm flux density of CMl from the publicly available CORNISH imag^ ( Pur cell et al. 



20081 ). and find a 6 cm flux d ensity of 10.7±0.4 mJy. We flt a free-free emission model 



as described in iHunter et al.l ( 120081 ) to the 6 cm-7 mm data (Fig. [3]), obtaining best-fit 
parameters Ne=9.5xl0^ cm-^ Te=10400 K, and size 0'.'4 (0.007 pc~1370 AU). This size 
agrees well with the deconvolved source sizes from two-dimensional Gaussian fits to our 3.6 
and 1 .3 cm images (~1230xll00 AU and ~1750xl060 AU, respectively). iBrogan et al. 
( I2OIII ) find an electron temperature of 7900 K and a density of 1.3x10^ cm~^ from analyzing 
H63a and H64a recombination line emission towards CMl; the lower values likely represent a 
contribution from a more extended, less dense component in the lower resolution EVLA data 
(beam 3'.'7 x 3'.'0, fitted source size 1'.'3). While its size would place CMl i n the HC HII regime , 



the electron density and modest recombi nation line width (17.6 km s ^ . iBrogan et al. 



are characteristic of UC HII regions (e.g. lKurtzll2005t iHoare et al.ll2007l : iMurphy et al 



2OIID 
2OI0I ). 



The ionizing photon fiux derived from our 3.6 and 1.3 cm observations (~1.5 xlO 
consistent with a single ionizing star of spectral type B1.5V. 



46 g-n 



IS 



While the deconvolved sizes of CMl and CM2 are quite similar (CM2: ~1750x820 AU), 
their spectra are very different. As shown in Figure |31 CM2 has a rising spectrum between 
3.6 cm and 7 mm (it is undetected in the 6 cm CORNISH image, 4cr limit 1.48 mJy). The 
spectral index between 3.6 and 1.3 cm is 0.67±0.14, derived from a Monte Carlo calculation 
of the spectral index that accounts for the uncertainties in the flux densities. The measured 
7 mm flux density (7.5±2.2 mJy, from the tapered and convolved image described above) is 
higher than predicted by extr apolation with g.^ fj i .sr-m (Fig. |3]). Th is 7 mm ex c ess ra ay be 
due to dust emission (see also IBrogan et al. II2OIII ). Alternatively, iKeto et al.l (|2008[ ) have 
suggested that HC HII SEDs with similar features can be flt purely as ionized emission 
from HII regions with steep density gradients. Mode ls of clumpy nebulae can al so reproduce 
intermediate spectral indices for HC HII regions (e.g. llgnace fc Churchwellll2004j ). In general, 



^Available at 



http : / / www. ast .leeds .ac.uk/ Cornish / public /releasel .php 
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an intermediate spectral index indicates that a range of opacities is present within a source. 
The deconvolved size of CM2 is significantly larger than the sizes derived assuming optically 
thick emission (Table E]). (Sub)arcsecond-resolution (sub)mm data are needed to constrain 
the relative contributions of dust and free-free emission in CM2. 

The natures of the other cm sources are only loosely constrained by our observations. 
The coincidence of C M3 with th e NE e xtension of the 24 yum emission suggests that it may 
be a MYSO; however, iDe Buizerl (120061 ) found that the 12 and 18 /xm MIR emission from the 
MYSO (and EGO) G35.2— 0.74 was dominated by emission from warm dust in the outflow 
cavity walls. The 24 /im ridge has approximately the same NE-SW axis as the bipolar 4.5 
/im lobes, so the possibility that the unsaturated 24 ^um emission NE and SW of CM2 traces 
an outflow cavity cannot be discounted. In this case, the ~symmetric spacing of CM3 and 
CM4 relative to CM2 might be suggestive of knots in an ionized jet. In any case, at least 
three (and possibly more) MYSOs are clustered at < 20000 AU scales, suggesti ve of a proto 
Trapezium such as seen in an increasing numbe r of m m (proto) clus t ers (e .g. iHunter et al. 
20061 : ICvganowski et allboO?! : iRodon et allboosl . Cll). ISewilo et all ^M) also found that 
all the HC HII regions in their sample were in pairs or multiples with other cm-A sources, 
though generally with wider separations (> 0.2 pc). 



3.2.10. G49.27-O.34 



Two cm continuum sources are detected coincident with the EGO, one resolved source 
that is detected at both 3.6 and 1.3 cm, and an unresolved source detected only at 3.6 
cm. The morphology of the resolved source, EGO G49.27— 0.34-CMl, is complex (Fig. [T]), 
with a bright, compact "head" (~2'.'2 x 1'.'2 ~0.06 x 0.03 pc), a knotty ring-like structure 
(~5'.'l X 4'.'6 ~0.14 X 0.12 pc), and more extended diffuse emission (dimensions measured 
from the 3.6 cm image). The 3.6 and 1.3 cm flux densities are consistent with optically 
thin free-free emission; the ionizing photon flux (~2.2xl0^^ s^^) corre sponds to a single 
ionizing star of spectral type BOV, in agreement with the estimate of iMehringerl (jl994( ) 
based on an unresolved 20 cm detection (~14" resolution). The spectral index limit for 
EGO G49.27— 0.34-CM2 (< 0.2) is also consistent with optically thin free-free emission, but 
the calculated ionizing photon flux is about two orders of magnitude lower (1.79x10^^ s~^)- 

Both cm sources are coincident with MIPS 24 /im em ission ; notably, neither is coinci- 
dent with the 44 GHz Cla ss I CH3OH masers detected by IC09l . This is of interest because 
De Buizer fc Vaccal ( 120101 ) obtained Gemini L and M band spectra for this EGO, and de- 
tected only continuum emission (no H2 or CO). Their slit crossed the bright multiband IRAC 
and MIPS 24 /im source near CM2, and passed near the northern 44 GHz maser (Fig.2 of 



- 14 - 



De Buizer fc Vaccal 120101) . The l ack of H2 emission in the Gemini spectrum and the 20 



cm detection by iMehringerl (119941 ) made this EGO a candidate for C lass I CH3OH masers 
excited by shocks from expanding HII regions ( iVoronkov et al.ll2010l ). but the large spatial 
offset complicates this interpretation. An outflow-perhaps driven by CM2, or by a MYSO 
undetected at cm wavelengths-may in fact be present, and responsible for the shocks excit- 
ing the maser emission. In this picture, shocked H2 or CO in the outflow would contribute 
to the 4.5 j»m emis s ion t o the south, near the Class I masers, but was either outside the 



De Buizer fc Vaccal (120101 ) slit or too weak to be detected in the Gemini spectrum. 



3.2.11. G49.42+O.33 

Two 3.6 cm sources are detected in the field, one (F G49.42+0.33-CM1) coincident with 
a multiband IRAC source and 24 ^m peak and onl y ~0 ?3 (> 4100 AU) from the eastern 6.7 



GHz CH3OH maser (G49.417+0.324) detected by lC09|. The other source, F G49.42+0.33- 
CM2, is extended and poorly imaged in our high-resolution VLA data; it is coincident with 
24 fim emission and an 8 //m nebula. If they are at the same distance as the EGO, the 
cm sources are > 0.7 pc (F G49.42+0.33-CM1) and > 1.8 pc (F G49.42+0.33-CM2) from 
the EGO position. It is unclear whether the cm sources are part of the same molecular 
cloud/star- forming complex as the extended 4.5 fim source. This EGO is also by far the 
most distant in our sample, so the cm nondetection at the EGO position corresponds to 
limits for the ionizing photon flux about an order of magnitude higher than for the other 
targets (Table |5]). 



4. Discussion 
4.1. Detection Rate: Comparison with Other Studies 

As stated in §21 we detect continuum emission at either 3.6 or 1.3 cm f rom 6 /14 (~43%) 



of the EGOs in our VLA survey sample. Including the four EGOs from the IC09I sample with 
deep continuum observations in the literature, all nondetections, brings the detection rate to 
6/18 (~33%)0 Since this sample is predominantly a 6.7 GHz CHsOH-maser-selected EGO 



^G23.96-0.11, G37.48-0.10, and G39.10+0.49 are undetected at 3.6 cm in VLA A-configuration observa - 
tions: la rms ~50 /xJy beam~^ (~60% higher than our median), Bgyn ~0'.'35 x 0'.'25 ( Bartkiewicz et allboool) 



G23. 01— 0.41 is undetected at 1.3 cm: 4tT hmit ~0.5 mJy beam ^(half our median), Ogyn (ICodella et al 



19971 ). We downloaded and reduced archival 3.6 cm VLA data (B configuration, project code AH868) for 
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subsample, the most relevant comparisons are searches for cm-A continuum emission towards 
Class II CH3OH masers. Two such recent st udies have sens i tivitie s and detection limits 
roughly comparable to our deep VLA survey. Pandian et al. (j2010) surveyed 20 6.7 GHz 



CH3OH masers selected from a blind single-dish maser survey, and iBartkiewicz et al.l (|2009[ ) 
targeted 30 6. 7 GHz CH.sOH rnasers sele cted from both blind and /■fi746'-targeted single dish 
surveys. The iPandian et al.l (I2OIOI ) and IBartkiewicz et al.l (120091 ) sampl es include 2 and 4 
EGOs from the ICOSi c atalo g, respectively, and 3 of the EGOs observed by IBartkiewicz et al. 
(I2OO9I ) are also in the IC09l sample. 



Pandian et al.l (120101 ) obtained shallow 3.6 cm and much deeper 1.3 cm VLA obser- 



vations (both with 9syn ~1")! a t 1.3 cm, they detect 30% (6/20) of the targeted masers 



in cm continuum emission. The iPandian et al.l ( 120101 ) data and our VLA 1.3 cm observa- 



tions provide comparable limits on central source properties; while their 1.3 cm observations 
are roughly twice as deep (la rms 0.12 mJy beam~^), their targets are significantly more 



distan t (median D=10 kpc, compared to 3.8 kpc for our EGO sample). IBartkiewicz et al. 



(120091 ) detect 3.6 cm continuum emission towards ~13% (4/30, with one additional ques- 
tionable associatio n) of their 6.7 GHz CH3OH maser sample, the distance range of which 
20 111 ) is similar to that of our target EGOs. Notably, in three sources 



(Bartkiewicz et a 




20111) 


Bartkiewicz et al. 


( 


2OO9I) 



emission with the continuum peak offset < 0'.'2 from the 6.7 GHz maser, similar to the 
continuum detections towards several EGOs (Fig. [H §3.2p . 

Considering our 1.3 and 3.6 cm detections separately for straightforward comparison, 
we detect ~31% of EGOs at 3.6 cm (4/13, 4/16 ~25% including literature data) and ~15% 
of EGOs at 1.3 cm (2/13, 2/14 ~14% including literature data). (These numbers exclude 
G49.27— 0.34, the only EGO in our sample without associated 6.7 GHz CH3OH masers.) In 
sum, the detection rat e for 1.3 cm continuu m emission towards EGOs from the |C09| sample 
is lower than that of IPandian et al.l (120 10[ ) towards 6.7 GHz CH3OH masers, though the 
limiting sensitivities in physical terms (e.g. ionizing photons s~^) are similar. In contrast, 
th e detection rate fo r 3.6 c m continuum emission towards these EGOs is higher than that 
of IBartkiewicz et al.l (|2009[ ) for their 6.7 GHz CH3OH maser sample. Considering the small 
number statistics for the detections, however, the detection rates for all three studies are 
roughly comparable. 



the remaining lC09l EGO. G10.34— 0.14. Unfortunately, imaging is compromised by the poor uv coverage of 
these data, and the rms is not constraining. 
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4.2. Nature of Cm Continuum Emission 

A number of diagnostic plots for distinguishing ionized emission from UC HII re- 
gions, HC HII region s, and MYSO winds and jets have been proposed (IHoare et al.l 120071 : 
Hoare &: Francdl2007l ). One key discriminant is the recombination line width, but such mea- 
surements are not available for most of our target EGOs. We can place our sources in Lj,(8 
GHz) versus Lhoi^iA yumj^ parameter space, albeit roughly, by using MIPS 24 /zm fluxes as a 
proxy for bolometric luminosity. 24 /im is the wavelength closest to the SED peak for which 
we currently have images with suflicient angular resolution to avoid source confusion. 



Mottram et al.l (120111 ) have explored using the MSX 21 fim broadband flux as a proxy 



for bolometric luminosity for YSOs and HII regions in the Red MSX Source (RMS) sample. 
For 613 sources with bolometric luminosities determined by fltting well-sampled SEDs, they 
find a mean Fftoi/FA/5'X2i^tm=21.9±1.9. While the ratio does not depend on source type 
(YSO V. HII region), the scatter is more than an order of magnitude (~7-200). For the 
MIPS 24 ^m filter, the same analysis of the RMS sample gives F;,o;/FM/p524Aim=29.5±1.9 
(J. Mottram, priv. comm.), with comparable scatter (~8-150). There is, however, some 
evidence that this mean ratio may underestimate the bolometric luminosities of EGOs. 

G19.01— 0.03 is a rare example of an EGO in which the "central" sou rce is a GLIMPSE 
point source, clearly resolved from extended emission in the IRAC images. ICllI observed this 
EGO at ~10,000 AU resolution in 1.3 mm continuum emission with the SMA, and detect 
only a single compact core, suggesting the observed MIR-mm emission arises from a single, 
embedded MYSO with a bo lometric luminosity of ~10^ Lq derived from SED modeling 
(iRobitaille et al.l l2006l . 120071 ). Extrapolating from the 24 /im fiux density using the mean 
Fbol/'PMIPS24,^lm ratio from the RMS data gives a bolometric luminosity of only 1.6x10^ 
Lq. It is reasonable to suppose that the SEDs of EGOs might peak at longer wavelengths 
than the MIR-bright population of YSOs and HII regions included in the RMS sample. 
If EGOs represent a younger stage, then the volume of surrounding dust that the central 
(proto)star has been able to heat will be much smaller than for an HII region, and so the 
average temperature will be lower. In this case, one would expect the RMS Fl,ol/FMIPS24:^im 
ratio to underpredict the bolometric luminosities of EGOs. For G19.01— 0.03, a range in 
Fboi/FMiPS2'inm of 8 to 150 corresponds to calculated bolometric luminosities of ~4xlO^- 
8x10'^ Lq; while still low, the upper end of this range approaches the bolometric luminosity 
derived from the SED fitting. 



^The quantity plot ted is the total bolometric luminosity, referred to by iHoare et al.l (|2007l ) and 
Hoare fc Francol (|20Q7I ) as L(IR). Wc estimate bolometric luminosities from 24 /im flux densities, and so 
use the term L;,o/(24 /im). 
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Considering the uncertainties in luminosity estimates and in the EGO 24 fim flux densi- 
ties, we a dopt t he following approach: For every EGO in our VLA sample with a 24 fim flux 
density in IcOqI . we calculate Lbo;(24^m) = S'24^m * 10"^^ *Au*x* 4 7rD^ * 2.48 x 10^^ where 
Lbo;(24 /im) has units of Lq, S24^m is the 24 //m flu x density fro m ICOOl in Jy, Au =2.8521 
xlO^^ Hz is the bandwidth of the MIPS 24 fim fllter (jCohenll2009l ). x is 8 or 150 (correspond- 
ing to the range in Fboi/FMiPS24:iim discussed above), and D is the distance to the source in 
kpc. Similarly, Lj,(8 GHz) is deflned a s 47rI )^S.':i.fif.m . These values are plotted in Figure H] 
along with data fromlHoare fc Francol ( 120071 ). In addition, data for low-mass sources from 
Table 5 of lAnglada et al. f ll998h are plotted for comparison Each EGO is plotted twice, 
with the upper estimate for L{,o/(24 fim) plotted as a flUed triangle, and the lower estimate 
for Lf,oi(24 fim) plotted as an open triangle for the same value of L,y{8 GHz). For EGOs 
with multiple cm counterparts coincident with MIPS 24 fim emission (G35.03+0.35 and 
G49.27— 0.34), all are plotted in Figure S] (sets of triangles at the same Lfeo/(24 fim) but 
different Lt/(8 GHz)). For nondetections, L,y(8 GHz) is calculated for the 4cr 3.6 cm upper 
limit. We consider Gil. 92— 0.61 a 3.6 cm nondetection and plot the Aa upper limit, rather 
than extrapolating from the 1.3 cm detection assuming a spectral index. All discussions 
of luminosities of course depend on distance estimates. For our EGO sample, we adopt 
the kinematic distances from C09. which are based on the new Galactic rotation model of 



Reid et al.l fl2009f ). 



Several salient points emerge from an examination of Figure |H (1) None of the EGOs 
extend into the upper right corner of the plot populated by UC HII regions. In particular, 
G35. 03-1-0. 35-CMl, for which we have good constraints on the electron density that place 
it flrmly within the UC HII category (§3X9]), has L^(8 GHz)~10^3 W Hz'^ towards the 
lo wer end of th e HC H II range. G35.03-f0. 35-CMl is less of an outlier in the line width plots 
of iHoare et al.l (120071 Fi g. 6). In line width/ size parameter space, the small recombination 
line width (17.6 km s'^ lBrogan et al. Iboili l and size (~0.007 pc) of G35.03+0.35-CM1 are 
consistent with the region occupi ed by UC HI I s, tho ugh G35. 03+0. 35-CMl has a narrower 
line width than any source in the iHoare et al.l ( 120071 ) sample. (2) EGOs detected at 3.6 cm 
generally lie in, or slightly above, the region of L,y{8 GHz) versus Lboi(24 ^m) parameter 
space occupied by winds and ionized jets. The EGO detecti ons to some extent flU in the 
gap between HC HIIs and wind/jet sources in the original iHoare fc Francol (120071 ) plot. 
However, as noted above, G35.03+0.35-CM1 and G49.27-0.34-CM1 fall in the HC HII region 
of Figure HI and their classiflcation as HCs is questionable. (3) The EGO nondetection 
upper limits all fall roughly in the region occupied by win d/jet sources, several orders of 
magnitude in L^{8 GHz) below the iHoare fc Francol ( 120071 ) UC and HC HII regions and 
more than an order of magnitude below the UC HII G35. 03+0. 35-CMl. (4) Even with 
highly uncertain, and possibly low, estimates for the bolometric luminosities, the properties 
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of the EGO population as a whole are consistent with MYSQ wi n d/iet sources and HC HIIs, 
and distinct from the low-mass sources from the lAnglada et al.l (119981 ) sample. 



In sum, our data definitively show that the vast majority (12/14^86%) of EGOs in our 
VLA sample are not associated with UC HII regions. The two EGOs that are associated 
with UC/C HII regions (G35.03+0.35 and G49.27— 0.34) also have other cm-A continuum 
counterparts, which may be at earlier evolutionary stages. The other EGO cm counterparts 
have spectral index limits consistent with either HC HII regions or winds/jets; their weak 3.6 
cm emission, relative to bolometric luminosity, is similar to well-studied wind and jet sources. 
In three EGOs (G24.94+0.07, G28.83-0.25, and G35.03+0.35), faint 3.6 cm emission is 
detected coincident with 6.7 GHz Class II CH3OH masers, indicating ionized gas very near 
the central MYSO. One important exception to the general trend of intermediate spectral- 
index emission is G11.92— 0.61-CMl, which may be an optically thick HC HII region. 



4.3. The Evolutionary State of EGOs 

As noted above, the majority of EGOs in our VLA survey (8/14~57%) are not associ- 
ated with detectable cm-A continuum emission in our deep observations. In evaluating the 
significance of cm-A nondetections towards MYSOs, there is always a degeneracy between age 
and mass: is the (proto)star not emitting ionizing radiation because it is young, or because 
it is not, and will never be, massive/luminous enough to do so? This degeneracy is to an 
important extent physical, and not simply an observational limitation: as a MYSO grows by 
accretion (whether competitively or otherwise), it will go through a stage at which its mass 
and luminosity are moderate, rapid accretion is ongoing, and its final mass is indeterminant, 
though limited by the mass reservoir of its natal clump. 

There is increasing evidence, from both modeling and observations, that accretion 
from a large-scale reservoir is a crucial element of the massive star formation process, 
and a key difference between clustered massive and isolated low-mass star formation (e.g . 



Galvan-Madrid et al.ll2009l : IPeters et al.ll201ll : IWang et al.ll2010l : ICalvan-Madrid et al.ll2010 



and references therein). So one way of observationally addressing the final mass of an actively 
accreting (proto)star, albeit indirectly, is to look at the mass of the clump-scale reservoir 
available for ongoing accretion. All but one of the EGOs in ou r VLA sample fall w ithin 



the area of the 1.1 mm Bolocam Galactic Plane Survey (B GPS; lAguirre et al. 



are associated with 1.1 mm sources in the BGPS catalog (iRosolowsky et al. 



2011 



2OIOI 



and 
^1 We 



^For the one source that fahs outside the BGPS survey area, Gi l. 92— 0.61, w e calcu late gas masses as 
described in the text but using the 850 /im flux density of 12 Jy from lWalsh et all |2003l) . R=100, and k^= 
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calculate gas masses from the 1.1 mm dust continuum emission 



M, 



gas 



4.79 X W-^'^RS^{Jy)D\kpc) 



(3) 



where is the flux density from the BGPS catalog co rrected by the recommended fac- 
tor of 1.5±0.15 ( lAguirre et al .1120111 : iDunham et al.ll2010l ). D is the distance to the source, 
B{u,Tdust) is the Planck function, R is the gas-to-dust mass ratio, and n^, is the dust mass 



Dunham et al. 



opacity coe fficient in units of cm^ g~^. We follow recent BGPS studies (e.g. 
2010l . I2OIII ) in adopting K27igh2/R=0.0114 cm^ g~^ (assuming a gas-to-dust mass ratio of 
100). For most of our sample, no information about the clump-scale dust or gas temperature 
is ayailable. We adopt two limiting dust temperatures, 16 and 28 K, bracketing the range 
(22±6 K) from lDunham et al.l (120111 ) for "Group 3" star-forming BGPS sources (those iden- 
tified with either an EGO or young RMS source). Clump masses calculated using equation [3] 
and the integrated fiux density of the associated BGPS source are listed in Tables [5] and |6l 
For three clumps that correspond to ATLASGAL 87 ^tm sources, the A TLASGAL mass 
(derived using estimates from NH3 observations; ISchuUer et al.l |2009| ) falls within our 
range of estimated clump gas masses. 

We do not attempt to correct the clump mass estimates for free-free contributions to 
the 1.1 mm fiux density for several reasons: (1) Many BGPS sources are large, and either 
encompass area outside the VLA 1.3 cm primary beam and/or extended cm sources for 
which fiux densities cannot be measured from our data. Hence we cannot correct for free- 
free contributions in a uniform way for all sources in our sample. (2) The uncertainties in the 
clump mass estimates are dominated by the uncertainty in the correction factor for BGPS 
catalog fiuxes (±15%) and in the assumed dust temperature. For the sources in our sample 
with the strongest free-free emission, G28.28— 0.36 and G49.27— 0.34, the total contribution 
of all cm sources within the VLA 1.3 cm primary beam corresponds to changes in the mass 
estimates of ~7% and <1%, respectively, assuming the steepest spectral indices allowed by 
our data (AM/M ~130/1840 Mq and ~45/5050 M0 for Td„«t=16 K, within the uncertainties 
quoted in Table [H]) . 

Figure [5^ shows histograms of masses estimated using Trf„st=16 K and the integrated 
fiux density from the BGPS catalog for EGOs with/without cm-A counterparts. For com- 
parison. Figure shows masses estimated using the 80" aperture fiuxes from the BGPS 
catalog. This provides a comparison of the mass at a uniform angular scale across our EGO 
sample, though the apertures are centered on the 1.1 mm peaks, and may be offset from 
the EGO positions. No trend in clump masses for EGOs with/without cm-A counterparts 



2.2 (see alsolCll|). 
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is seen in Figure [5^ or b. Adopting Tdust='28 K for all sources does not change the relative 
positions of cm detections/nondetections, only shifts the absolute mass scale. If we instead 
assume that sources with cm-A continuum detections are warmer, the mass range of clumps 
associated with EGO cm detections and nondetections still largely overlaps (Fig. [St,d). 



For BGPS clumps associated with HII regions, iJohnston et al.l (120091 ) found a correla- 
tion between the clump mass and that of the ionizing stars. This indicates a relationship 
between the mass of a clump, and the mass of the massive stars that it forms. This idea, 
combined with the importance of accretion from large-scale reservoirs to massive star forma- 
tion, suggests that the EGOs in our survey sample have the potential to continue accreting, 
and eventually become massive stars that will ionize their surroundings (i.e. form UC HII 
regions). As shown in Figure there is no apparent trend in clump masses for EGOs 
with cm-A detections as compared to cm-A nondetections, as might be expected if the pres- 
ence/absence of cm continuum emission was related to reservoir-and so to final stellar-mass. 
In this picture, the paucity of cm-A continuum detections towards EGOs is a measure of their 
youth. 

Additional support for this explanation is derived from the SED fitting of G19.01— 0.03, 
a cm nondetection. This MYSO is actively driving a bipolar outflow, traced by Cl ass I 



CH3OH masers, SiO emission, and high-velocity, coUimated ^^CO and HCO"*" lobes (jC09 



ICl j ). An active outflow implies active accretion, which is expected to alter the properties 
of the central (proto)star, affecting its output of ionizing radiation. Based on the SED 
modeling, the central (proto)star of G19.01— 0.03 h as a luminosity of ~10^ Lq, mass of > 10 



Mq, and is accreting at a rate of ~10 ^ M© yr ^ (ICllI ). Figure [6^1 shows the stellar radius 
and temperature for the best-fit models 
of the IRAC, MIPS, and SMA data from 



Xbest datapoint <3) to the SED composed 
Gill . Both swollen, very low temperature (~5000 



K) and hotter (>10^ K) models are allowed; however, even the highest stellar temperature 
in Figure ^ falls below the 30 , 000 K threshold for "significant" Lyman continuum output 
suggested by iHoare fc Francol (120071 ) . Figure [6]3 shows that if the lo w angular resolution 



single dish) (sub)mm flux density estimates at 870 /im and 1.1 mm (jSchuUer et al.l 12009 



Rosolowsky et al.ll2010[ ) are substituted for the SMA data, all the well-fit models are pushed 
to large radii and low temperatures^. 



Longmore et al.l (120091 ) report the detection of spatially extended, optically thin free- 



free emission in the fields of six 6.7 GHz C H.sOH masers previously reported as 3.6 cm 
continuum nondetections (IWalsh et al.l Il998l ). and caution that the evolutionary state of 



^°For both plots, SEDs we re fit assuming a min imum uncertainty for all flux densities of 10% to mirror 
recent single-dish studies (e.g. lPandian et al.ll201Cll ). which differs somewhat from the errors assumed bv lCll 
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CH3OH maser sources may be misconstrued if extended free-free emission is missed in high- 
resolution interferometric observations. We beheve that our cm continuum study of EGOs 
is robust against such misclassification for several reasons. (1) The sen sitivity and inaagiri g 
fidelity of our VLA observations are better than the ATCA survey of IWalsh et al.l ( 119981 ): 
though we cannot recover fluxes for sources more extended than ~10", our data are sensitive 
to structures up to ~20" in angular scale This is illu strated by ou r 3.6 c m detection of 
the HII region F G28.83-0.25-CM1, which was missed bv IWalsh et"aD (Il998h (2) We have 
examined GLIMPSE and MIPSGAL 24 /xm images for the iLongmore et al.l (120091 ) targets, 
and the extended HII regions reported manifest themselves as structured, multiband IRAC 
sources and/or 8 fim n ebulae associated with bright 24 fim emission (see also discussion of 
"diffuse red clumps" by lBattersby et al.ll2010l ). This is consistent with the MIR morphologies 
of our "field" cm sources (Fig. H]), including those overresolved in our observations (and so 
candidates for being po ssibly "extended a nd ov erlooked" ) . (3) In several cases, the extended 
HII regions reported by lLongmore et al.l (120091 ) are offset from the CH3OH maser positions, 
and the masers are more closely associated with HC HII regions or with 3 mm emission 
attributed to dust. So, as iLongmore et al.l (120091 ) note, the presence of more evolved HII 
regions in the vicinity of 6.7 GHz CH3OH masers may simply be due to the clustered nature 
of massive star formation. In fact, we see several potential examples of this in our data. 
In these cases, the high-resolution GLIMPSE and MIPSGAL MIR images are helpful in 
identifying evolved sources and assessing their likely relation to (less evolved) sources of 
interest (gSD- 



5. Conclusions 



Our exceptionally deep VLA 3.6 cm survey (median rms ~30 /iJy beam~^) provides 
convincing evidence that most EGOs are young MYSOs, prior to the development of a 
UC HII region. No continuum emission is detected at 3.6 or 1.3 cm (median rms ~0.25 
mJy beam~^) towards the majority (8/14 ~57%) of EGOs in our survey sample. Only 
two sources (2/14 ~14%) are associated with optically thin cm-A emission from UC/C HII 
regions (one of which is the only one of our targets without an associated 6.7 GHz CH3OH 
maser). Each of these EGOs exhibits cm-A multiplicity, with morphological evidence that 
a less-evolved cm continuum source may be the driving source of the outflow traced by 4.5 
fim and Class I CH3OH maser emission. With the exception of the UC HII regions, most of 
the cm-A EGO counterparts are weak (< 1 mJy), unresolved at the ~1" resolution of our 
survey, undetected at 1.3 cm, and characterized by intermediate spectral indices consistent 
with either HC HII regions or ionized winds or jets. Three of these weak, unresolved 3.6 cm 
sources are coincident with 6.7 GHz CH3OH maser emission, indicative of ionized gas very 
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near the central (proto)star. Likely to be an optically thick HC HIT region, one EGO cm 
counterpart is detected only at 1.3 cm and is associated with a hot core and H2O and 6.7 
GHz CH3OH maser emission. The predominance of very weak cm-A continuum detections 
in our survey indicates the promise of the greatly enhanced continuum sensitivity of the 
Expanded Very Large Array (EVLA) for understanding the early stages of massive star 
formation, including the earliest appearance of ionized gas around accreting MYSOs. 

In sum, the cm wavelength properties of EGOs are indicative of an early stage of massive 
star formation, when photoionizing feedback from the central MYSO is not yet a significant 
influence on the (proto) cluster e nvironment . The presence of 4.5 fim emission. Class I 



CH3OH masers, and SiO emission (IC09l : ICllI ) indicates that these EGOs are actively driving 
outflows, and so (presumably) actively accreting. Notably, two of the EGOs with weak 
or no cm emission in our survey (Gil. 92— 0.61 and G19.01— 0.03) have coUimated, high- 
velocity molecular outflows, in contrast to some recent theoretical predictions that ionizing 



radia tion will dominate outflow dynamics in massive star-forming regions (e.g. iPeters et al. 



201 ll ). Already sufficiently massive to excite 6.7 GHz CH3OH maser emission, the surveyed 
EGOs are associated with significant clump-scale gas reservoirs. These clumps can provide 
sufficient material for sustained, rapid accretion as the EGOs continue to grow from MYSOs 
into massive stars. 
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Table 1. Observational Parameters 



Source Name 


J2000.0 Coordinates'' 


Phase Cal. 


Syn. Beam'' 


3.6 cm 


1.3 cm 


uv 




„ (h m s) 


5 {° ' ") 


(1.3 cm) 


"x "[°] 


rms 


rms 


range? 












(mjy beam"-*^) 


(mJy beam^^) 


(kA) 


GlO.29-0.13 


18 08 49.3 


-20 05 57.3 


1820-254 


1.45x0.86(0.38] 


0.136 


0.32 


>10 (3.6cm) 


Gil. 92-0.61 


18 13 58.1 


-18 54 16.7 


1820-254 


1.45x0.85[1.51] 


0.052 


0.25 




G18.67+0.03 


18 24 53.7 


-12 39 20.0 


1832-105 


1.25x0.91[-2.03] 


0.031 


0.23 




G18.89-0.47 


18 27 07.9 


-12 41 35.5 


1832-105 


1.25x0.91[-3.77] 


0.037 


0.25 


>25 (3.6cm) 


G19.01-0.03 


18 25 44.8 


-12 22 45.8 


1832-105 


1.24x0.93[-3.72] 


0.029 


0.26 


>10 (3.6cm) 


G19.36-0.03 


18 26 25.8 


-12 03 56.9 


1832-105 


1.24x0.92[-4.07] 


0.029 


0.26 




G22.04+0.22 


18 30 34.7 


-9 34 47.0 


1832-105 


1.20x0.95[-0.24] 


0.031 


0.27 




G24.94+0.07 


18 36 31.5 


-07 04 16.0 


1832-105 


1.14x0.95[-1.39] 


0.036 


0.25 


>10 (1.3cm) 


G25.27-0.43 


18 38 56.9 


-07 00 48.0 


1832-105 


1.16x0.96[-1.88] 


0.030 


0.26 




G28.28-0.36 


18 44 13.2 


-04 18 04.0 


1851+005 


1.12x0.95[6.75] 


0.053 


0.24 




G28.83-0.25 


18 44 51.0 


-03 45 49.0 


1851+005 


1.09x0.95[7.52] 


0.029 


0.23 


>10 (3.6,1.3cm) 


G35.03+0.35 


18 54 01.0 


+02 01 20.0 


1851+005 


1.05x0.96[8.92] 


0.028 


0.20 




G49.27-0.34 


19 23 06.7 


+14 20 13.0 


1924+156 


0.96x0.89[-7.35] 


0.135 


0.18 




G49.42+0.33 


19 20 59.3 


+14 46 48.0 


1924+156 


0.93x0.88[-17.83] 


0.027 


0.17 





^Phase tracking centers of the VLA observations; in some cases, these may differ slightly from the lCOsI positions. 
''At 1.3 cm; the 3.6 cm data were convolved to this resolution, see J2] 



Table 2. EGO Detections: Observed Properties of Centimeter Continuum Sources 





Source Name 






3.6 cm 








1.3 cm 










J2000.0 Coordinates'' 


Peak Intensity'' 




Size'^ 


J2000.0 Coordinates'* 


Peak Intensity'' 




Size<^ 






„ (h m s) 


5 (° ' ") 


(mJy beam~^) 


(mJy) 


("x "[P. A. °]) 


„ m S-) S (° ' ") 


(mJy beam~^) 


(mJy) 


("x "[P. A. °]) 


EGO 


G11.92-0.61-CM1 






<0.21 






18 13 58.10 -18 54 20.5 


1.07(0.25) 






EGO 


G24.94+0.07-CM1 


18 36 31.572 


-07 04 16.93 


0.53(0.04) 








<1.01 






EGO 


G28.28-0.36-CM1 


18 44 13.33 


-04 18 04.3 


0.21(0.05) 








<0.96 






EGO 


G28.83-0.25-CM1 


18 44 50.737 


-03 45 49.26 


0.33(0.03) 








<0.92 








-CM2 


18 44 51.08 


-03 45 48.3 


0.14(0.03) 








<0.92 






EGO 


G35.03+0.35-CM1 


18 54 00.49098 


+02 01 18.292 


12.34(0.03) 


13.76(0.05) 


0.36x0.32(111] 


18 54 00.4902 +02 01 18.327 


10.30(0.20) 


12.07(0.38) 


0.51x0.31(157] 




-CM2 


18 54 00.6498 


+02 01 19.32 


1.29(0.03) 


1.49(0.05) 


0.51x0.24[21] 


18 54 00.647 +02 01 19.53 


2.88(0.19) 








-CM3 


18 54 00.766 


+02 01 22.82 


0.45(0.03) 








<0.81 








-CM4 


18 54 00.523 


+02 01 15.60 


0.23(0.03) 








<0.81 








-CMS 


18 54 00.45 


+02 01 21.2 


0.13(0.03) 








<0.81 






EGO 


G49.27-0.34-CMl<^ 


19 23 06.91 


+14 20 16.7 


2.44(0.14) 


73(9) 


9.5 


19 23 06.91 +14 20 16.7 


2.63(0.18) 


67(7) 


9.1 




-CM2 


19 23 06.61 


+14 20 12.0 


0.61(0.14) 








<0.71 







to 

''Tlie number of significant figures indicates tlie formal uncertainty in tiie centroid position obtained from a two-dimensional Gaussian fit, or (for sources <5 tr), a one-pixM'unccrtainty 
in the position of the peak (see also 



''All upper limits are 4(t. For detections, the quoted error is the formal error from the two dimensional Gaussian fit, or the Icr rms in the image, depending on whether the source was 
detected with sufficient signal-to- noise to fit (see Q. 

•^Measured from two-dimensional Gaussian fits; formal errors from the fits are given in parentheses. 

'^Source properties measured using polygonal aperture photometry program developed for the CORNISH survey. The quoted uncertainty in the peak intensity is the Icr rms in the 
image. The quoted size is the equivalent source diameter, defined as the diameter of a disk with the same solid angle as the source polygon. See also [[3] 
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Table 3. EGO Nondetections: Observed Centimeter Limits 



Source Name 


3.6 cm 


1.3 cm 




Peak Intensity* Peak Intensity'' 




(mjy beam ^) 


(mJy beam ^) 


EGO GlO.29-0.13 


<0.54'' 


<1.27 


EGO G18.67+0.03 


<0.12 


<0.94 


EGO G18. 89-0.47 


<0.15 


<0.98 


EGO G19.01-0.03 


<0.12 


<1.04 


EGO G19.36-0.03 


<0.12 


<1.03 


EGO G22.04+0.22 


<0.12 


<1.06 


EGO G25.27-0.43 


<0.12 


< 1.03 


EGO G49.42+0.33 


<0.11 


<0.68 



"■All upper limits are 4(T. 

noise level is elevated due to image artifacts 
caused by copious complex, extended emission in the 
field (and beyond the half-power point of the primary 
beam) that is poorly imaged in our observations. 



Table 4. Observed Properties of Additional Centimeter Continuum Sources in EGO Fields 



Source Name 






3.6 cm 






1.3 cm 






Ang. Sep.'^ 




J2000.0 Coordinates'' 


Peak Intensity'' 




Size 


J2000.0 Coordinates'* 


Peak Intensity'^ 


" s^ 


Size 


from EGO 




a m =) 


(5 (° ' ") 


(mjy beam"-^) 


(mJy) 




^ (h m s) 5 (o / //) 


(mJy beam~^) 


(mJy) 




(") 


F G10.29-0.13-CM1 


18 08 50.6 


-20 06 00 


c 


C 


>21" 


18 08 50.6 -20 06 00 


c 


e 


>21" 


>io 


F G18.67+0.03-CM1 


18 24 52.6056 


-12 39 20.005 


4.29 (0.03) 


6.35(0.07) 


0.83 X 0.63 [41]f 


18 24 52.607 -12 39 20.02 


3.7(0.4) 






16 


F G19.36-0.03-CM1 


18 26 25.628 


-12 03 49.43 


0.42(0.03) 


0.99(0.09) 


1.39xl.06[55]f 




<1.03 






8 


-CM2 


18 26 26.379 


-12 04 19.90 


0.48(0.03) 








<1.16 






24 


F G25.27-0.43-CM1 


18 38 59.6635 


-07 00 55.15 


1.02(0.03) 








< 1.23 






40 


-CM2 


18 38 57.576 


-07 01 14.78 


0.31(0.03) 








< 1.16 






28 


F G28.28~0.36-CM1S 


18 44 15.11 


-04 17 56.1 


77.32(0.05) 


576(3) 


7.7" 


18 44 15.11 -04 17 56.2 


69.6(0.3) 


510(11) 


8.9" 


30 


-CM2 


18 44 09.802 


-04 18 00.12 


1.51(0.07) 


2.0(0.1) 


0.69x0.54[42]f 




<1.47 






51 


-CM3 


18 44 16.498 


-04 18 08.86 


1.03(0.06) 


3.0(0.2) 


1.53xl.33[176]'' 




<1.69 






50 


F G28.83-0.25-CM1 


18 44 51.3 


-03 45 22 


C 


c 


~27" 










>13 


F G49.42+0.33-CM1 


19 20 59.81 


14 46 49.4 


0.14(0.03) 








<0.68 






11 


-CM2 


19 20 56.4 


14 46 44 


C 


c 


~20" 










^30 



00 

''Measured as described in [|3]for different source types. For sources fit with Gaussians, the number of significant figures indicates the formal uncertainty in the centroid piteition 
from the Gaussian fit. 



''For detections, the quoted error is the formal error from the two dimensional Gaussian fit, or the la rms in the image (depending on whether the source was detected with 
sufficient signal-to-noise to fit, ijSj, except as otherwise noted. 

■^AU upper limits are 4(t, where tr is the rms at the location of the source in the primary-beam-corrected image (for sources far from the pointing center, this will differ from 
the rms at the EGO position). 

■^Angular separation between the centimeter emission and the EGO position from lCo"^ (minimum angular separation for extended cm sources). 

^Extended source poorly imaged by the interferometer, so no fluxes are reported. 

'Deconvolved source size ("x "[P. A. °]) as determined by fitting a single 2D Gaussian component. 

^Source properties measured using polygonal aperture photometry program developed for the CORNISH survey. The quoted uncertainty in the peak intensity is the Icr rms 
in the image. The quoted size is the equivalent source diameter, defined as the diameter of a disk with the same solid angle as the source polygon. See also ^ 
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Table 5. EGO Cm Nondetections: Derived Physical Properties and Limits 



Source Name 


Distance'' 
(kpc) 


(mpc) 


3.6 cm'' 
(AU) 


^LyJ 

(s-1) 
xlO** 


(mpc) 


1.3 cm'' 

!<= 
(AU) 


(s-1) 


Clump 
Td=16K 
(Mq) 


mass"^ 
Td=28K 
(Me) 


L^(8 GHz) 
(W Hz-i) 


L6<,,(24 fimf 
min max 
(L©) (L©) 


EGO GlO.29-0.13 


2.19 


<0.39 


<80 


<2.5 


<0.22 


<46 


<0.7 


12201^20 




<3.1xl0" 






EGO G18.67+0.03 


4.98 


<0.42 


<86 


<2.8 


<0.44 


<90 


<2.5 


1660ll« 




<3.6xl0" 


6.9x10^ 


1.3x10" 


EGO G18.89-0.47 


4.49 


<0.42 


<86 


<2.9 


<0.40 


<83 


<2.1 


5350lg7Q 


2520t«o 


<3.6xl0" 


7.2x10' 


1.3x10^ 


EGO G19.01-0.03 


4.20 


<0.35 


<72 


<2.0 


<0.39 


<80 


<2.0 


11401^^0 


^40+110 


<2.5xl0" 


4.4x10^ 


8.2x10^ 


EGO G19.36-0.03 


2.43 


<0.20 


<42 


<0.7 


<0.22 


<46 


<0.7 


109011* 




<8.5xl0'0 






EGO G22.04+0.22 


3.62 


<0.30 


<62 


<1.5 


<0.34 


<70 


<1.5 


1650l^~ 


780ll« 


<1.9xl0" 


3.1x10^ 


5.8x10^ 


EGO G25.27-0.43 


3.86 


<0.32 


<66 


<1.7 


<0.36 


<73 


<1.6 




4501^00 


<2.1xl0" 


2.1x10' 


4.0x10^ 


EGO G49.42+0.33 


12.29 


<0.98 


<202 


<16 


<0.92 


<189 


<11 


344Q+900 


1620t«o 


<2.0xl0'2 







'"From lcpgl . 

''All quantities are upper limits corresponding to a 4a nondetection in our VLA images. 

Clump masses calculated from the int egrated flux density measurement in the BGPS catalog l lRosolowskv et aI]|201Cll ). multipled by 
a factor of 1.5±0.15 as recommended by lAguirre et all 1I2OIII) . The quoted uncertainties include the uncertainty in this multiplicative 
factor. See ^431 

''Estimated from the 24 /im flux density from lcoil see i|4.2l ■ ■ ■ indicates no 24 fim flux density in i COSj . either because a clear 
counterpart is not present, or because confusion/blending precludes measuring a reliable flux for the EGO counterpart. 

''Twice the radius calculated in equation [2] ([|3TT}- Assumes optically thick emission. 

'Assumes optically thin emission, i|3.1l 



Table 6. EGO Cm Detections: Derived Physical Properties and Limits 



Source Name 



Distance^ 
(kpc) 



3.6 cm'^ 



1.3 cm " 
1' N 



Lyc 



Clump mass'^ U(8 GHz) Lj,o,(24 ^tm)<= 
Td=16K Ttj=28K (W Hz^l) min max 



(mpc) (AU) (s-i) (mpc) (AU) (s-i) (Mq) (Mq) 



(L©) (L©) 



xlO* 



EGO G11.92-0.61-CM1 3.80 >1.7 <0.42 

EGO G24.94+0.07-CM1 2.99 <0.7 0.52 

EGO G28.28-0.36-CM1 3.29 <1.6 0.36 

EGO G28.83-0.25-CM1 5.03 <1.1 0.70 

-CM2 5.03 <1.9 0.45 

EGO G35.03+0.35-CM1 3.43 -0.1 ■ • • ' 

-CM2 3.43 0.7 1.01 

-CMS 3.43 <0.6 0.55 

-CM4 3.43 <1.3 0.40 

-CM5 3.43 <1.9 0.30 

EGO G49.27-0.34-CM1 5.55 -0.1 ■•• ' 

-CM2 5.55 <0.2 ■•• ' 



<86 <2.89 0.36 



73 



1.66 



1100'' 



108 4.51 <0.27 <56 <0.97 340^7° 



75 



93 



82 
61 



-1+600 



154 



15.3 2070^ 



208 16.7 0.53 109 3.65 207OI330 
114 5.04 <0.28 <58 <1.03 2070^330 

360 
330 



2.58 < 0.28 < 58 < 1.03 2070; 



1.46 <0.28 <58 <1.03 2070Z 



2140 
17.9 



222 5050 



+900 



<2.35 5050 



+900 
-830 



490'' 



870lli° 
1580^280 
15801^80 

980i}^° 

980li^° 
2380t«» 
2380l«° 



<3.6xl0" 
5.7x10" 
2.7x10" 
1.0x10^2 
4.2x10" 
1.9x10^^ 
2.1x10^2 

6.3x10" 

3.2x10" 

1.8x10" 

2.7x10^"' 
2.2x10^2 



1.1x10-^ 2.0x10" 
8.9x10^ 1.7x10^ 



6.5x10^ 1.2x10'* 

1.5x10^ 2.8x10"* 

1.5x10^ 2.8x10* 

1.5x10^ 2.8x10"* 

1.5x10^ 2.8x10"* 

1.5x10=* 2.8x10"* 

2.4x10^ 4.5x10"* 

2.4x10^ 4.5x10"* 



''From lCOgl . 

'^Spectral index a, oc i/", calculated between 3.6 and 1.3 cm from our VLA data. 
'^AU limits correspond to 4(7 upper limits in our VLA images. 

'^Clump masses calculated from the integrated fl ux density measurement in the BGPS catalog llRosolowskv et al]|201ol) . multiplcd by a factor of 



00 
00 



1.5±0.15 as recommended bv lAeuirre et all ll201lh . The quoted uncertainties include the uncertainty in this multiplicative factor. See ^4.31 



"Estimated from the 24 fim flux density from [C09l . see i|4.2l ■ ■ ■ indicates no 24 fim flux density in lC09l . either because a clear counterpart is 
not present, or because confusion/blending precludes measuring a reliable flux for the EGO counterpart. 

'Twice the radius calculated in equation[2l ( i|3.1l l. Assumes optically thick emission. 

^Assumes optically thin emission, i|3.1l 

'^Gll.92— 0.61 is the on ly EGO in our sam ple that falls outside the BGPS coverage. The masses listed here were calculated from the 850 ^m 
SCUBA flux reported by 'Walsh ct al.' (20o3) (12 Jy), assuming Tdust ^ ^ov the BGPS sources and K850pm=2.2 (interpolated from the values 
tabulated by Ossenkopf & Hcnning 1994). 

'Not calculated, since the radius estimate assumes unresolved optically thick emission ( i|3.1l l. and the spectral index indicates that the emission 
is optically thin. 
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Fig. 1.- 
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(f) G49.27-0.34 





Fig. 1.— Cm-A EGO detections. Three-color Spitzer GLIMPSE images (RGB: 8.0, 
4.5, 3.6 fim) overlaid with contours of VLA cm continuum emission (blue) and MIPS 24 
fim emission (yellow). All panels are centered on the VLA pointing center except (e), 
which is centered on the position of G35 .03+ 0.35-CM2. Positions of 6.7 GHz Class II 
and 44 GHz Class I CH3OH masers from IC09l are plotted as red diamonds and magenta 
crosses, respectively. The VLA beam is shown at lower left in each panel (Table [1]). 
The VLA images shown have not been corrected for the primary beam, but all mea- 
surements were made from corrected images, (a) G1L92— 0.61: VLA 1.3 cm continuum 
contours at (3,4) xo" = 2.54E-4 Jy beam~^ and CARMA 1.4 mm continuum (black con- 
tours), levels (4,5,10,20) xa = 4.3 mJy beam~^ (t he C ARMA beam is shown at lower 
left). The H2O maser from iHofner fc Churchwelll (119961 ) is plotted as a green triangle. 
MIPS 24 fxm contour levels: 600,900,1200,1500,1800,2100 MJy sr~i (b) G24.94+0.07: VLA 
3.6 cm continuum contour levels (4,7,12) xo" = 3.57E-5 Jy beam~^. MIPS 24 /im con- 
tour levels: 300,600,900 MJy sr^^ (c) G28.28-0.36: VLA 3.6 cm continuum contour lev- 
els (3,4) xa =5.25E-5 Jy beam'^ MIPS 24 //m contour levels: 900,1200,1500,1800 MJy 
sr-i (d) G28.83-0.25: VLA 3.6 cm continuum contour levels (3.5,4.5,7,9) xa = 2.92E-5 
Jy beam-i. MIPS 24 fim contour levels: 600,900,1200,1500,1800 (e) G35.03+0.35: VLA 3.6 
cm contour levels (4,7,12,36,108,324) xa = 2.78E-5 Jy beam'^ MIPS 24 fim contour levels: 
900,1200,1800 MJy sr'^ (f) G49.27-0.34: VLA 3.6 cm contour levels (4,7,10,13,16) xa = 
1.35E-4 Jy beam-^ MIPS 24 fim contour levels:600,900,1200, 1500,1800 MJy sr^^. 
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Fig. 2. — Field Cm-A Continuum Sources. Three-color Spitzer GLIMPSE images (RGB: 
8.0, 4.5, 3.6 //m) overlaid with contours of VLA 3.6 cm continuum emission (blue/cyan), 
BGPS 1.1 mm emission (white), and MIPS 24 emission (yellow). All panels are cen- 
tered on the VLA pointing center. Positions of 6.7 GHz Class II and 44 GHz Class I CH3OH 
masers from lC09l are plotted as diamonds and magenta crosses, respectively. The VLA beam 
is shown at lower left in each panel (Table [T]). The VLA images shown here have not been 
corrected for the primary beam, but all measurements were made from corrected images, 
(a) GlO.29-0.13: VLA 3.6 cm continuum contour levels: (4,7,10) xa = 1.36E-4 Jy beam'^; 
BGPS 1.1 mm contour levels: (10,40,80) xcr = 0.02 Jy beam"^; MIPS 24 fim contour levels: 
900,1200,1500,1800,2100,2400 MJy sr'^ (b) G18.67+0.03: VLA 3.6 cm continuum contour 
levels: (4,40,120) xa = 3.06E-5 Jy beam"^; BGPS 1.1 mm contour levels: (10,20) xcr = 
0.03 Jy beam"^ MIPS 24 contour levels: 600,1200,1800 MJy sr'^ (c) G19.36-0.03: 
VLA 3.6 cm continuum contour levels: (4,12) xcr =2.93E-5 Jy beam^^; BGPS 1.1 mm con- 
tour levels: (10,20,40) xcr = 0.03 Jy beam^^; MIPS 24 fim contour levels: 600,1300,2000 
MJy sr-i (d) G25.27-0.43: VLA 3.6 cm continuum contour levels: (4,7,21) xa = 2.95E- 
5Jy beam-i; BGPS 1.1 mm contour levels: (10,20) xa = 0.018 Jy beam'^; MIPS 24 /im 
contour levels: 75,150,600 MJy sr~^ (e) G28.28— 0.36: VLA 3.6 cm continuum contour levels: 
(4,10,20,50,150,250,450,650,850,1250) xa = 5.25E-5 Jy beam~i; BGPS 1.1 mm contour lev- 
els: (10,20,50,80) xa = 0.015 Jy beam"^; MIPS 24 ^m contour levels : 900,1200,1500,1800 
MJy sr-i (f) G28.83-0.25: VLA 3.6 cm continuum contour levels: (4,5,7) x a =2.86E-5 
Jy beam-i; BGPS 1.1 mm contour levels: (10,20,50,100) xa = 0.015 Jy beam^^; MIPS 24 
fim contour levels: 600,1200,1800 MJy sr^^ (g) G49.42+0.33: VLA 3.6 cm continuum con- 
tour levels: (4,5,7) xa = 2.69E-5 Jy beam"^; BGPS 1.1 mm contour levels: (5,10,15) xa = 
0.025 Jy beam-i; MIPS 24 fiia contour levels: 150,300,450 MJy sr'^ (h) Same as (g); only 
the 4 and 5 a VLA 3.6 cm contours are visible within the s maller field of view. The BGPS 
rms is a function of Galactic longitude ( lAguirre et al.ll201ll Fig. 11). 
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Fig. 3.— SEDs for the cm sources EGO G35.03+0.35-CM1 (left) and G35.03+0.35-CM2 
(right). The dot-dashed hne in the left panel is the fit described in §3.2.9[ In the right panel, 
lines of a =0.67 (dot-dashed) and of ai la =0.14 (dotted), extrapolated from the 3.6 cm 
flux density, are shown ( §3.2.91 ) 
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Fig. 4. — Plots of Lj,(8 GHz) v. bolometric luminosity Lfeo/(24 /iin) for EGOs detected 
at 3.6 cm (t op) and for EGO nond etections (bottom). UC HIIs, HC HIIs, wind, and jet 
sources from iHoare fc Francol (120071 ) are plo tted as X's, so l id dot s, open circles, and open 
squares, respectively; low-mass sources from lAnglada et al.l (Il998[ ) are plotted as stars. For 
each EGO cm counterpart detected, the datapoint corresponding to our lower estimate for 
Lboi(24 /im) is plotted as a right-facing open triangle, and the datapoint corresponding to 
our upper estimate for Lfeo/(24 //m) is plotted as a left-facing solid triangle, with the two 
connected with a line (see §4.2p . Nondetections are plotted as downward-pointing triangles, 
with a solid triangle for the upper and an open triangle for the lower L^oii^A yum) estimate, 
joined by a line. 
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Fig. 5. — Top: Histograms of the masses of clumps associated with EGOs calculated from 
the integrated flux dens ities (a) and 80" aperture flux densities (b) in the BGPS catalog 
(iRosolowsky et al.l |2010| ) for a dust temperature of 16 K. For clarity, only nominal mass 
values are plotted. Histograms of all clump masses (EGO cm detections and nondetections) 
are plotted as sohd lines, clumps associated with EGO cm detections are plotted as hatched 
histograms. Bottom: Same as top, except Tdust='^8 K is used for EGOs with cm detec- 
tions. Cm detections and nondetections are plotted as diagonally and horizontally hatched 
histograms, respectively. 
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Fig. 6. — Stellar radius and temp erature for mode l s well-fit (y ^ — xlest P^^ datapoint <3) to 



the G19.01-0.03 SED using the iRobitaille et all (120061 120071 ) model grid and fitter (black 
dots; the greyscale shows all models in the grid), (a) Fits to the datapoints used in ICll 



(IRAC, MIPS 24 and 70 ^m, and SMA 1.3 mm); (b) fits to the same MIR datapoints but 
with ATLASGAL 870 fim and BGPS 1.1 mm data in place of the SMA 1.3 mm point. 



